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ABSTRACT 


Fatigue~crack growth data for 2024-~T3 aluminum was used 
to evaluate the accuracy of Forman's equation. Forman's 
equation was found to predict from one-third to twice the 
observed cycles to a given crack length when applied to 
multiple cyclic stress ratios. An alternative approach 
was developed, based on the observed data, which is 
applicable to analog or digital fatigue-crack growth 
prediction. Data acquisition based on constant increment 
of crack length has been found to yield inadequate definition 
of the crack~growth curve in the low-cycle range for the 
successful application of the method. An expression was 
developed which extends the application of the crack-closure 
concept to the entire range of cyclic stress ratio. An 
effective root-mean-square stress was defined which shows 
indications of good correlation with the magnitude of the 


crack~growth curves. 
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I. INTRODUCTION 


Modern aircraft, particularly military aircraft, must be 
designed for high performance and long service life. The 
attempt to achieve higher performance in terms of increased 
payload, range, speed, or maneuverability implies a reduc- 
Bec FrOSSsaircrait welght, or in theystructural fraction 
of that weight. This constraint results in structural 
commonents which carry higher loads per unit area and places 
greater demands on their structural integrity. 

The design values for the safe-life of non-inspectable 
components and the inspection intervals for other components 
are predicated upon the service life and the predicted load 
spectrum specified by the contracting agency. The aircraft 
components are sized to meet these design values using 
fatigue prediction techniques. The estimates of remaining 
life and maximum safe inspection intervals of operational 
aircraft are frequently re-evaluated in light of the opera- 
tional load spectrum using similar methods. A factor of 
Saeeuy, Or factor of uncertainty, is used in the’se calcula— 
tions by which the expected life of a given structural 
configuration is reduced. This factor of safety (FS) 
directly reflects the accuracy of the prediction method used. 
Large FS, indicating ae accuracy, results’ in ekcessive 


structural weight when compared to the optimum and penalizes 


Chew ceratit performance throughout the service life. 


i. 





Optimistic estimates of the accuracy of a method lead 
ultimately to costly retrofits and reduced aircraft 
availability. ‘The long term economic efficiency of an 
aircraft is therefore dependent upon the accuracy of the 
fatigue prediction methods available to the designer and 
operator; thus considerable motivation exists for the 
improvement of fatigue life prediction accuracy. 

The fail-safe design philosophy, which admits the 
existence of a crack in the structure, is one of the fatigxe 
prediction approaches in current use. Crack growth predic 
CLONS eine units Of cycles from an initial size to its lengin 
when it precipitates catastrophic failure, are an essential 
element in the fatigue analysis of today's airplanes. 
Inspection periods and rework intervals can be prescribed 
from this information. This thesis reports an investigatin 
of data and theory used widely, if not exclusively, in 
making crack growth predictions for aluminum alloys. 
Specifically, Hudson's fatigue crack growth data [1] for 
2024-T3 aluminum alloy was examined in order to establish 
confidence limits on Forman's equation, which appears to b:. 
the most widely accepted crack growth prediction procedure 
used. Finally, an alternate method for expressing crack 
growth behavior is proposed to improve the accuracy of cry 


length predictions. 


I 





II. BACKGROUND 


The fatigue life of aircraft components is difficult 
to predict due to the many factors influencing the structural 
response to unsteady loading. These factors may be grouped 
into four categories; material characteristics, environmental 
effects, structural configuration, and service loading. 

The fatigue response to unsteady ioads is statistical 
in nature and is dependent upon the time history of the 
loading and material response, where material response is 
defined to be cycles to a given crack length. The response 
history may be divided into three stages. The first stage, 
termed crack initiation, is characterized by sub-microscopic 
cracks. This stage has been shown to be sensitive to a 
variety of microscopic damage mechanisms. The intermediate 
stage, termed crack propagation, is characterized by tne 
relatively slow growth of a macroscopic crack (a few grain 
diameters) under the influence of the ordered external load. 
Final failure occupies a small fraction of the fatigue life 
of the material and a very high crack growth rate exists 


during this stage. 


A. MINER'S RULE OF CRACK INITIATION 

The obvious complexity of the problem has led investigators 
to make certain simplifying assumptions and/or approximations. 
Miner [2] proposed a linear damage accumulation over the crack 


initiation stage of the fatigue life. This proposal defines 


15) 





damage as the fraction of the initiation life in cycles, Ny; 


oh load level. An increment of damage 


would then take place with each cycle at the i°" 10a4a level 


expended at the i 


according to 1/N, or for n, cycles, the damage would be 
n,/N, . Since the accumulation of damage was assumed to be 
linear, Miner's Rule predicts the fatigue life of a 


component as 


L = : ny (1) 
and damage as 
yap (2) 
a i 
At the completion of the crack initiation stage, 
RL, ee (3) 


a 


Fatigue life data is presently available as So Vs Ne 
curves where So 1s the alternating component of stress and 
Np is the number of cycles to complete failure of the 
specimen rather than formation of a first crack. Application 
of Miner's rule to random load tests have yielded damage 
predictions at specimen failure considerably different from 
unity. Some refinements to Miner's approach to damage 
Pecdictlomsiave sparen wae form of weighting functions on each 


m 
term or variable powers on each term such as D = & (n,/N, ) : 
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These modifications as well as competing damage hypotheses 
have ofttimes resulted in improved prediction accuracy for 

a given set of data; however, the Miner summation of damage 
continues to be the most widely accepted due to its simplicity 
and the minimal improvement in accuracy offered by its more 


complicated alternatives. 


B. CRACK PROPAGATION 

After the crack is formed, fail-safe design practice of 
crack-growth predictions requires a knowledge of the crack 
growth rate and the load under which the cracked structure 
would fail, both of which are not available from any damage 
theory. Many investigators, as summarized in Ref. 3, have 
attempted to correlate the crack growth per cycle to the 
gross section stress and the crack length, leading to 


propagaticn models of the forn, 


da/aN = f(S,a,A) (4) 


where 2a is the total crack length, 5S is the gross section 
stress and X is a correction factor which compensates for 
finite width, rivet forces, etc. Paris and Erdogan [4] 

concluded that the stress intensity factor (SIF), which is 
a function of the same variables, was the parameter which 


best correlated with the crack growth rate observed. 


ats, 





C. STRESS INTENSITY FACTOR 

Mathematically, the stress intensity factors are undeter- 
mined coefficients arising from the solution of the field 
equations for any cracked geometry. The normal stress at 
the crack tip, under moderate static loading, is given by 


this solution as 


K K K 
Gey = —e £,,(8) + Eg, ,(0) + fh, , (8) (5) 
My 1G IG 


where Ky > K,> and K. are the opening mode, in-plane shear 
mode, and out-of-plane shear mode SIF respectively, r is 
the radial coordinate referenced to the crack tip and 98 is 
the angular coordinate measured counter-clockwise from the 
advancing crack tip. The stress intensity factors may 
therefore be interpreted as parameters that reflect the 
PeGieit | i enmot ne sliressetilecidsinga body resulting from 
the introduction of a crack. The magnitude of the stress 
intensity factor is dependent upon the geometry of the body 
Penvainiing Gres crack, the size,and location of. the crack, 
and distribution and magnitude of the external loads on the 
Dedvwee lic con! Wcurablon mCSterecadily treated 1s the sheev 
specimen containing a central crack of length 2a under 


Mode 1 loading, for which 


ae Si fra | (6) 
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where dw is a width correction factor. The criterion for 
failure in the presence of a crack-like defect > cha sc racks 
growth to failure (instability) will occur whenever the 
applied stresses, S, as given by equation (6) exceed some 
eritical condition specified by K, = Geo Paris and Ergodan 
[4] concluded from the preceding argument that this condition 
may be extended to fatigue and applied to the slow growth 


phase. 


D. FORMAN'S EQUATION 
Forman [5] proposed a model for crack propagation which 


utilized the range of the SIF as suggested by Paris [6], 
AK = K -K. (7) 


to account for the cyclic loading typical of the fatigue 
environment. This model included the instability condition 
for failure as a singularity such that the growth rate 


approaches infinity as Ko approaches Kor’ Forman's 


AX 


equation is written as 


da _ CCAR 
aN CI=-R)K,), - Ak (8) 


where C and n were assumed to be undetermined material 
COnstvaelivcmand R is thetraitilo of minimum stress to maximum 


stress (loading ratio). 


aby 





E. HUDSON'S DATA 

Hudson [1] conducted a most extensive study of the 
material response of two widely used aircraft alloys, 
2024-T3 and 7075-T6 aluminum, to evaluate Forman's equation 
when applied to a large number of maximum stress conditions 
and loading ratio combinations. (The data obtained in 
Hudson's investigation was used as the data base for this 
study.) 

The specimens used in Hudson's study were nominally 
12 inches in width, 35 inches in length, and 0.090 inches in 
thickness. A notch 0.10 inches in length by 0.01 inches in 
width was cut into the center of each specimen by an elec- 
trical discharge machining technique. A reference grid of 
0.050 inch spacing was photographically printed on the 
surface of each specimen to mark intervals in the path of 
ene sea cks 

"Patigue-crack growth was observed through 10-power 
microscopes while illuminating the specimen with stroboscopic 
light. The number of cycles required to propagate the crack 
to each grid line was recorded so that crack-propagation 
rates could be determined. Approximately two-thirds of the 
crack-propagation tests were conducted to failure. The 
remaining one-third were stopped before failure, and the 


elac Coomapecincis Were Uscdein residual—-static strength tests." 


INational Aeronautics and Space Administration Technical 


Note D-5390, Effect of Stress Ratio on Fatigue-Crack Growth 
in 7075-T6 and 2024-T3 Aluminum Alloy Specimens, by C. 


Michael Hudson, p. 5, August 1969. 
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The results of Hudson's fatigue-crack-growth tests on 
the 2024-T3 specimens are presented in Table I (reproduced 
from Ref. 1). Hudson obtained the fatigue-crack propagation 
rates graphically from the crack-growth curves defined in 
Table I. The crack-growth curves reproduced from Ref. 1 
and shown as Figure 1 are typical. The results were 
presented as plots of the common logarithms of the experi- 
mental value of the crack-propagation rate as functions of 
the range of the SIF. Three models were then fitted to the 
data by least-squares techniques and compared to the data. 
hucures. = (rom met. ©. Is ty¥pilceal of tne mevnod used for the 
presentation of data and the evaluation of the quality of 
fit of empirical equations. Hudson concluded that Forman's 
equation produced an excellent fit for both materials whereas 
the other equations evaluated "...produced good correlation 
with the test data except at the high growth rates for the 
7075-T6 alloy" on the basis of graphical representations of 
which Figure 2 is typical. This conclusion at first appears 
to be justified on the basis of plots such as Figure 2, 

Howe VeiIe = Gie, Lapleue analyse ls concerned With tre accuracy 
with which a crack-propagation model predicts the number of 
cycles required for the crack to attain a given length. This 


is the basis upon which the fit should be evaluated. 


ing 





a ee Ob VELOEMENT 


A. EVALUATION OF FORMAN'S EQUATION 
Hudson used the tangent form of the width correction 


PveL@l’s 


W Ta (9) 


in evaluating the material constants, C and n, in Forman's 
equation. The differential equation resulting from the use 
of this form of width correction in Forman's equation may 
be integrated (Appendix A) using the small angle approximation, 
tan® ~@ . This approximation is accurate to within three 
percent for crack lengths of less than 0.8 inches for 
Hudson's specimens. 

The predicted interval of cycles between crack lengths, 
AN, was calculated using Hudson's fit of Forman's equation 


Pp 


and the integrated form of Forman's equation becomes 


Da O- 
AN. = 2¢ |. 2 (her : Smax Vo, _ pat (rer _ Snax V2 | 
p i O n-2 n-3 n- n-3 
(10) 
where ~~. @1_R) 
L n 
CE(1-R)S, oo] 
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and AND is the predicted number of cycles required to 
propagate at half-crack of length ao to lengthias The ratio 
of the predicted to the observed cycle increment is presented 
Piet toe Ie Wwoemecm cross “disirepanctesmaxist for almost tevery 
case. 

The load cases were also fitted individually in order to 
determine the cause of the poor predictions reported in 


Table II. The expression for AK used in this analysis was 


= + ae 
218 ( Saspte Seat “Aw Vf Fi (11) 


where Brown and Srawley's [7] finite width correction for a 


centrally cracked sheet, 


s eI, 2ay\e 2a,3 
was chosen on the basis of its reported accuracy. The test 
data for each load case were represented by a power series 
of the form 


C5 a 
+ poN + ... + p_N (1B) 


which yields a very good approximation of the observed data 
as indicated by Figure 3. An analytic crack-growth rate 
expression was obtained by performing a differentiation of 


equation (13) with respect to cycles leading to 


eal 





e,-—l e,-1 
aay _ 1 2 
an > P e4N + Poe5N 


e 
+e ot Pom a (14) 


Equations (11) and (13) were applied to each load case and 
a least-squares fit of Forman's equation to the data was 
performed. Figure 4 is typical of the individual fits thus 
obtained. | 
| A regular, periodic variation of the calculated crack- 
growth rate about the fitted line was observed which is 
believed to indicate that Forman's equation may be an 
incorrect or incomplete functional form. These periodic 
variations appear as scatter when Forman's equation is called 
to fit multiple load cases. An apparent lower bound of 
crack-growth rate may also be observed in Figure 4. The 
failure of Forman's equation to represent this threshold 
levels probably accounts for the largest portion of prediction 
error. Just as for the integration performed for Hudson's 
fit, consistent prediction errors were noted for the fits 
of the individual load cases but of smaller magnitude. The 
Perlogue, Variat@on Ofsthefdatasabeuu thePfitted line Vand the 
apparent lower threshold of crack-growth rate led to the 
development of a new material response model. Equation (13) 
was not considered tobe a practical alternative to Forman's 
equation in view of the explicit appearance of the variable 
N, which is unknown in the practical application of any 
prediction model. Series reversion was likewise discounted 


due to non-integer exponents resulting from the fit. 


Pah (24 





B. MATERIAL RESPONSE FUNCTION 

A material response function was developed in which the 
response was assumed to be a piecewise continuous function 
of time, loading conditions, and crack geometry. The 
function was also assumed to possess piecewise continuous 
time derivatives. This assumption was based on the satis- 
factory representation of the data achieved by the power 
series oe ation ef the fatigue data, equation (13), 
which has an infinite number of continuous time derivatives. 

Previous investigators have assumed that the fatigue- 


crack propagation rate, da/dN, has the functional form 
da/aN = f(S(N) ,a(N) ,A(N) ) (5) 


where the variables have the meanings previously defined and 
are functions of time or, alternatively, cycles. Althcugh 
the relationship between the independent variables N and t 


is linear, as expressed by 
N = ot (16) 


where w is the circular frequency of loading, a distinction 
in the growth rates da/dt and da/dN exists. The crack-growth 
rate da/dt is an instantaneous rate which must be integrated 
over the period of a cycle in order to obtain the time 


average rate, da/dN. 


23 





Paris [6] concluded that the crack growth rate was best 
represented as a function of the range SIF as given by 
equation (11). The desired function must predict crack 
growth rate under constant amplitude and random loading with 
equal precision since constant amplitude loading is a special 
case Of random loading. The stress range is therefore a 
constant for either load case taken on a cycle-by-cycle basis. 
This argument applies equally well to any other choice of 
the characteristic stress to be used in calculating the SIF. 
The expression for crack-growth rate could then be represented 


equally well by the functional form 


aA = £,(a(N),A(N)) (17) 
OT 
Ga = £ (x(N)) (18) 
aN O 
where 


x = d, (asW)\/a(N) 


since A is a function of crack geometry only. 
Introduction of the requirement for higher order deri- 
vatives can be satisfied by a linear summation of the assumed 


form 


dx 


a’ x 
= ¢) an + .,, + —_— (19) 


re 
| 


Q_xX 
“2 me 
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where these, ="c,(S). Equation (19) isea linear differential 
equation in N with constant coefficients and thus has a 
sOlution of exponential functions. Paris [6] noted that the 
erack growth rate is a monotonically increasing function 

of the 43-8! intensity factor. The erack lerneth, a, 
inereases monotonically with time and the width correction, 
X, increases monotonically with crack length, therefore the 
rate of change of the variable, x, with respect to cycles 
(time) must increase monotonically with x. This type of 
response is typical of positive feedback systems, which 
implies that the roots of the characteristic equation 
obtained from equation (19) are real. A mathematical model 
of the desired material response function is therefore the 


Prony series 


x = : Dos 4 explb,.N] (20) 
Which satisfies all requirements heretofore set forth. 
This form of the material response function may be recast 


to the differential 


aX*x _ 1 ax dex , ae ca J 
2 ee Sk 


——=— = — [ec ste 
k 1 AN aN 
(al) 


aN Wi 


when the®*functional form of the Cy is determined. This 
Giprereny won eduacsonm would find wtility if programmed for 


analog or hybrid computation using non-linear elements 
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eliminating the excessive digital computer time requirements 
typical of random loading fatigue calculations; however, it 

is beyond the scope of this thesis. This approach to fatigue 
prediction calculations would also allow the use of discretized 
flight load data without intermediate processing. 

Thaemdatay sampling frequency used by Hudson, based on 
equal increments of crack length, provided insufficient data 
in the low cycle range of the specimen lives. The crack 
growth rate progresses from a nearly constant value at small 
crack lengths to much larger values within a small increment 
of crack length requiring a more frequent data sampling 
frequency in order to achieve a successful analysis of the 
exponential form of the fatigue response function through 


direct analysis. 


C. EFFECTIVE STRESS 

The extension of equation (20) to account for random load 
sequences requires that the constants be a function of the 
loading conditions. The representation of the piecewise 
continuous response function, x'(t), as the cycle average 
function x(N) leads to the representation of the piecewise 
continuous crack-growth rate, da(t)/dt, as the time average 
crack-growth rate, da(N)/dN. Consistency requires that the 
loading also be represented as a time average variable. 
Paris [16] suggested the utility of the Root-Mean-—Square 


(rms) gross section stress, as defined as 
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T 9 
ip Sm Ge che. (22) 
0 


~” 
HI 
Hh 


rms 


where S(t) = Sm + Sa sin wt 


and Sm is the mean gross section stress during the cycle, 
Sa is the amplitude of the gross section varying stress 
during the cycle, w is the circular frequency and T is the 
period of the cycle. 

The entire magnitude of the stress described above dces 
not contribute to crack propagation. Elber [8] has noted 
that plastic yielding at the crack tip results in a compres- 
Sive residual stress upon unloading. This compressive stress 
at the crack tip delays the opening of the crack, Elber 


defined the crack closure parameter, U, as 


meer OD 
max min 


where So is the gross section stress at which the crack 


p 
opens. He found that U seemed to be a linear function of 


loading ratio for 2024-T3 aluminum as represented by 
U= 0.4R + 0.5 (24) 
DICcomiemerackmprolaratvtom Oceurs While the crackeis 


closed, an effective rms gross section stress may be defined 


as 


Zi 


an 





ems - : 2. oT (25) 
where ce Sue Ta om 

e m a ak Op 
subject to 5. > 0 


The nature of the sinusoidal function is such that the limits 
of integration can be shifted to write (25) more conveniently 
in terms of the damage effective portion of the cycle as 


oor dv Co leage 2) Oe eeas so keads CoO 


yes gts 
Ny Peary q~ oT 
2 _ 2k 4 at. 
“rms. =f 6 7 iu 0 T (26) 
iD Ee iy ge 
ome Ot Ne rr. 


where f/T is the fraction of the cycle during which the 


crack is open and is defined by (See Appendix B) 


cos7+(1-20) 
TI 


ic 


f/T 2 
Elber's equation for U(R) leads to effective cycle 
fractions which are undefined for load ratios less than -1.25. 
This unnatural behavior in the neighborhood of a subset of 

Hudson's data motivated a search for a more satisfactory 
representation of the crack closure phenomenon. Figure 7 
Us Ba schiematterw or the foading eonditions for increasingly 
negative load ratios. Reasoning from this figure indicates 


that, in the limit, as 4 approaches zero (R approaches 


ax 
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negative infinity), the fraction of the period over which 
the crack is open approaches zero. This physical constraint 


is satisfied by an equation of the form 


a. 
7~5- (28) 


where A and B were assumed to be undetermined material 
constants. Substitution of equation (28) into the defini- 


iron get f/T ,eaquation (2/7), leads to 





U= 3% [1 - cosa (g4.)1 (29) 


Elber's data [8] was fitted to the function for the 
effective cycle fraction, equation (28), as shown in Figure 
6. The form satisfies the physical constraint .and yields 
a smaller Moa error than equation (24) when compared in 
Figure 7. Figure 7 is a comparison of the functional 
representation of U(R) as developed in this study with that 
of Elber. Both forms appear to be incomplete in that neither 
Porimcencl ders syhnesenrects "Ol mbsolure. Loadglével or crack 
geometry on the crack closure phenomenon. S&S was 


rmS g 
caleulated on the basis of equations (25) and (238). 


D. HYBERBOLIC REPRESENTATION OF FATIGUE-CRACK GROWTH DATA 
On the basis of the arguments presented with regard to 
the exponential form of the material response function and 


the observed similarity of the exponential and hyperbolic 


ae 





7 
- 


curve forms, Hudson's fatigue-crack growth data was fitted 
to the general form of the second order equation, 
ye 2 as 
x + BXN tOCN@e+ Dxat EN —= & + (30) 
Solving for N, Tablelll is a listing of the coefficients of 
the equation 


2 
N = d,x t dy + 4/d2x tdyx + de 5 (31) 


3 


ordered with decreasing 8 which resulted from the fit 


agate 
Of equation (30). 

Figure 8 indicates the quality typical of these coeffi- 
cients. The variation in sign which occurs in Table III 
reflects both the translation of the x axis (or the N 
ordinate of the center of the fitted hyperbola) and the 
frequent occurence of parabolic, rather than hyperbolic, 
coefficients. The parabolic coefficients (sign of A = sign 
of C in equation (30)) resulted from the paucity of data in 
the low-cycle range of the specimen fatigue life. The 
appearance of parabolic coefficients is unacceptable from 
ee Daysicalm@standpoine as Che parabolic form implies that 
cycles to a given crack length is not a single valued 
function, as shown schematically in Figure 9a. The trans- 
lation of the center of the assumed hyperbola was an 


unavoidable result of the presence of the original machined 


30 





notch. Figure 9b is a schematic representation of the 
effect of load level and the presence of the artificial 


notch on the translation of the x axis. 


E. NORMALIZATION OF THE DATA 

A unified fatigue prediction method may be formulated 
Beyecomnbinune the capabilities ofsthe damage accumulation 
method in estimating crack initiation with the fracture 
mechanics approach applied over the crack propagation and 
failure stages. This usage demands the normalization of 
all load cases to a common reference crack length which 
should be the minimum detectable crack length defining 
initiation. This normalization was accomplished for each 
(hyperbolic) load case by estimating the cyclic interval 
between Hudson's reference point, typically 0.10 inch, and 
the assumed initiation points, 0.003 inch. The estimate was 
made on the basis of equation (31) using the d, previously 
calculated and was then applied as an increment to the 
observed cycles to a given crack length for an unnotched 
specimen. The modified data sets were again fit to equation 
(30) yielding estimated coefficients, d,*, of the unnotched 
sheet which are hypothesized to correlate with ernse 
Figures 10 through 14 indicate the degree of this correlation 


and suggest that the appropriate functional relation may be 


S i (32) 


Bul 





An analytic expression for the crack-growth rate, da/dN, 


may be obtained by the chain differentiation 


ANI as Lone (ax) (au (33) 
da ax 
Oxy aL A dA 
where nee —— ( 5 i 2 Ga? (34) 


1 

= (2a.xt+d,,) 

aati GN = SS (35) 
dx Fdyx tds 


This expression for da/dN approaches the constant (threshold) 
crack-growth rate observed in Figure 4 at smail crack lengths 
and satisfies the instability condition given that equation 
is). hnalsas zero. 


The zero of equation (34) yields x and leads to 


ins 


Kins . max ins (36) 
where the instability condition has been assumed to take 
place at the instant in the last cycle at which the load 
has attained its maximum value. Figure 15 depicts the 


variation of the ratio K ay © . This apparent reduction of 


a er 


tne level of the’stress intensity factor for unstable crack 


Bre 





propagation is attributed to dynamic effects on propagation 
and has the effect of minimizing the influence of the 
Singularity condition in Forman's equation. 

The data fits to equation (30) yielded no real solution 
to equation (34) for some load cases. These load cases 
appeared to be the dynamic tests which were terminated well 
before failure and therefore contained insufficient 
information regarding the shape of the response function in 


the high-cycle region to give the correct physical response. 


Bee) 





IV. CONCLUSIONS 


Hudson's [1] fatigue-crack growth data for 2024-73 
aluminum was used to evaluate the accuracy of Forman's 
equation. An alternative approach to the representation 
of the propagation process was developed and i Sit Ae An 
expression was developed which extends Elber's crack-closure 
concept to the possible load ratio range. The following 
conelusions can be drawn from this study. 

1. Forman's equation does not appear to satisfactorily 
predict crack propagation in 2024-T3 when evaluated on the 
basis of the observed data. 

2. Elber's crack-closure concept was successfully 
extended to the entire load ratio range in 2024-T3 by the 


relation 


U = sin (st) 


3. The development of a crack-growth rate expression 
based on empirical analysis of the observed data met with 
limited success due to the coarse data resolution in the 
low-cycle range. 

Wmnece Stress intemsiey factorsavp erack» instability 
under dynamic loading conditions was observed to decrease 
with decreasing values of effective rms gross-section 
stress, asymptotically approaching a value of approximately 


one-third the value measured under static loading. 
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IV. RECOMMENDATIONS FOR FUTURE RESEARCH 


The material response of 2024-T3 aluminum has been shown 
fo be capable of representation by a continuous response 
mone. LO ic accCUMmmMbatlOn Ofmamnumber of Gest sets with 
more frequent data sampling in the low-cycle phase of the 
crack propagation life would allow a determination of a 
material response model for which statistical reliability 
limits could be determined. <Any such model must be based 
upon the observed data, as compared to derived data, for 
the reliability limits to be meaningful. 

The development of the exponential response model will 
allow the determination of the governing differential equation 
for crack propagation under dynamic ioading by the application 
of parameter identification techniques. Such a development 
would allow the calculation of the expected crack length 
under random loading by the use of analog, rather than 


disicaly eompucatwone devi ces. 
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APPENDIX A. INTEGRATION OF FORMAN'S EQUATION 


Forman's equation is written 


da _ __C(AK)" 7 STS 
aN Git pre AV 


where the stress intensity factor, K, is defined as 


k 4 sav/ra (n=2) 


and the range of the stress intensity factor is written 
AK = AS\N/ra (A-3) 
The tangcens form so: the width correct#on factor, 


or = ae tan a (A-4) 
may be used, allowing the small angle approximation. This 
apeuemimaviOnw tan) Om——80., yields! errors of Jess ‘than three 
percent for ratios of crack half-length to width of less 
thgnea/15. 

Eamlitane the integration te 1th s rangem@allows AK Co be 


expressed as 
AK = AS\fa (hes) 
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The cyclic stress ratio, R, is defined 


5 | 
R = = (A-6) 
min 
Therefore AS = Sh) Sa (A-7) 


Forman's equation can be rewritten as 


n n _n/e 

= “ COl-R) Siax (A-8) 
Gok, -s07% 
ian ke ) 


for the limitations imposed. 
Treating this expression as a differential as proposed 


by Forman [5] 


(1-R)(K.. - Sal’) 
c(1-R)"s_. "a 


which may be integrated between crack half-lengths ao and a 


172 
Wo. eae Zz cep) axe oN 
P C(1-R)"S_ ZI ee 7 
C80) 
Catiev | m2 eoue the indicated anvegration yields 
e—-n } 3a a 
( ) K a 2 Sea te! 2 
n 
C(1-R) Say eS — a 





Simplification of equation (A-11) leads to 


o-n 
o Ror Snax 


p ih O n-2 n-3 


where 


(1-R) 


C(1-R)"S z 


ax 


i 
J 
NO 
= 
1 
La) 


(A-13) 


and AN, is the cyclic increment predicted to elapse during 


the extension of a crack of half-length a 


to a half-length 


of a when Subjected to constant amplitude loading conditions 


R and ae 


38 








APPENDEX BY THE DERIVATION OF THE RELATIONSHIP BETWEEN 
THE DAMAGE EFFECTIVE CYCLE FRACTION AND 
THE CRACK=~CLOSURE PARAMETER 
Crack propagation occurs only when the fatigue crack 
is open, allowing the definition of a damage effective 
portion of a cycle (period). This fraction of the period 
is bounded by the limits t,/T and to/T representing the 
instants of crack opening and crack closing, respectively. 


The damage effective fraction is then 
f/T = t5/T - t,/T (B-1) 


Hesse Lilies are Symmevric about the meximum load instant 
for constant amplitude loading. Assuming a loading function 


of the form 


S(t/T) = Sin ou Se S ie 2 iu (B-2) 


the limits of the damage effective fraction of the cycle 


may be written 


W/4 - f£/(2T) (B-3a) 


ays 


and 


to/T = W/4 + £/(27) (B-3b) 
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The crack opens during the loading portion of the period 
when the effective gross section stress is zero. The crack- 
closure concept modifies this stress by the inclusion of 
the residual stresses resulting from local tensile yielding 
at the crack tip. Defining Sop to be the gross section 


Stress at which the crack opens leads to 


S. (oy) = Sa = Sop + So Sin (2nt/T) (B-4) 


The instant of crack opening is therefore defined by 


Ss (t,/T) = BLS S58 ae et (21t,/T) =50 (B-5) 


oT 


S - $ 
(AE —*) = sin [2nr(1/4 - f/2T)1 (B-6) 


a 


Trigonometric reduction yields 
=o m 
(oA) = ges (n/t) - (B-7) 


AnGmunc andre weitectuuvemmractilon, Or the cycle is}defined 
in terms of the loading and closure parameters 


Le 
£/t & cost*/n (—OP_™) (B-8) 


S 
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Rewriting Equation (B-8) in terms of ae and 
yields 
3 5 
3 = ( max min ) 
“7 Op 2 
2 
which simplifies to 
(Seo) 
f-/2 = cos t/t ja __max OP 
S =a 5 
max min 
This becomes 
a -1 
f/T = cos “/tm [1 - 2U] 


with Elber's [8] definition of the crack-closure 


U, : 


> 
Op 
min 


S 
ma X 


max 


Solving Equation (B-12) for U yields 
i oe al 
= 5 = ces Tia 


Therefore, the relation between f/T and U is 


2 (aly 


sin ap 


4 


Miu 


(B-9) 


Os Sa1 00) 


(B-11) 


parameter, 


=p 


(B-13) 


(B=14) 





and 


i 
2 


e/T = 2/n sin? (U) (B-15) 
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TABLE II. Analysis of the prediction accuracy of Forman's 
equation for small cracks in 2024-T3 aluminum. 
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TABLE III. Coefficients, di, resulting from the fit of the 
fatigue-crack growth data to the general 
_ second order equation’ 
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Figurescem=buadcon'’s comparison of fatigue-crack mrowth rate 
with Forman's equation (From Ref. 1). 
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Figure 4. Comparison of fatipfue-crack frowth rate with 
Forman's equation for a single load case (typical). 
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